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Active  Feedback  Control  of  Combustor  Dynamics 
with  Time  Delay  and  Noise 

C.  Seywert  G.  Isella  F.E.C.  Culick 

California  Institute  of  Technology 
1200  East  California  Boulevard,  MS  205-45 
Pasadena,  CA  91125,  USA 


Abstract 

Longitudinal  pressure  oscillations  in  a combustion  chamber  are  studied  theoretically  by  means  of  a 
low  order  model  obtained  by  systematic  reduction  from  a complete  representation.  The  formulation  is 
based  on  the  derivation  of  a generalized  wave  equation  that  accommodates  the  effects  of  mean  flow, 
combustion,  noise  and  control  action.  By  using  spatial  averaging,  the  equations  describing  the  dynamics 
of  the  chamber  are  reduced  to  a set  of  coupled  ordinary  differential  equations.  The  form  of  the  resulting 
equations  is  particularly  convenient  for  model  reduction  and  for  introducing  feedback  control  terms,  while 
retaining  all  physical  processes. 

Active  feedback  control  of  longitudinal  instabilities  is  then  introduced  by  using  the  same  formulation, 
rewritten  in  state-space  form.  The  broad  objective  of  control  simulations  covered  here  is  to  investigate 
in  a unified  fashion  various  aspects  of  the  problem.  These  include  the  influences  of  noise,  parameter 
uncertainties,  unmodeled  modes  and  time-delay. 

A criterion  is  derived  that  guarantees  the  stability  of  the  controlled  closed-loop  system  in  the  presence 
of  those  quantities.  The  particular  controller  used  here  is  based  on  a standard  LQR  design,  but  any  design 
technique  can  be  used  as  long  as  the  stability  criterion  is  fulfilled. 


1 Introduction 


The  need  to  extend  the  operating  range  of  combustion  systems  due  to  more  stringent  performance  demands 
(lower  emission  levels,  reduced  vibration  tolerances,...)  has  led  to  increased  research  activity  on  control 
of  combustion  instabilities.  Traditionally  these  pressure  oscillations  have  been  dealt  with  through  passive 
control  techniques.  However  since  these  techniques  are  inadequate  under  varying  operating  conditions,  active 
control  methods  have  become  more  widespread.  McManus,  Poinsot  and  Candel  (1993)  gives  a review  of  the 
different  approaches  taken  nearly  a decade  ago.  Considerable  progress  has  been  made  since  that  time, 
particularly  with  investigations  of  the  suitability  of  known  control  strategies  applied  to  combustor  dynamics. 

Combustion  systems  are  highly  complex  and  are  naturally  represented  by  infinite-dimensional  models.  So  as 
to  reduce  the  order  of  the  controllers  several  researchers  have  made  use  of  an  approximate  finite-dimensional 
representation  of  the  flowfield  (Culick  1989)  as  the  basis  for  their  controller  design  (Hathout,  Annaswamy, 
Fleifil  and  Ghoniem  1998,  Dowling  1997,  Haddad,  Leonessa,  Corrado  and  Kapila  1997,  Koshigoe,  Komat- 
suzaki  and  Yang  1999,  Hong,  Yang  and  Ray  2000).  Since  any  model  is  prone  to  uncertainties  in  the  para- 
meters, the  controller  needs  to  be  robust  to  changes  in  those  parameters.  Haddad  et  al.  (1997)  achieves  this 
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robustness  using  a Lyapunov  function  framework  while  Koshigoe  et  al.  (1999)  achieves  robust  performance 
by  using  adaptive  control.  Other  authors  (Hong  et  al.  2000)  use  an  /Too  technique  to  reach  the  same  goal. 

The  typical  characteristics  of  combustors  make  designing  a controller  a challenging  task.  The  model  used  in 
the  design  is  normally  a considerably  reduced  representation  of  the  real  system;  this  causes  large  uncertainties 
in  the  parameters.  In  addition  to  that,  there  are  very  intense  internal  noise  sources  due  to  gasdynamics  and 
combustion  and  significant  time  delays  that  reduce  the  stability  margin.  Scaling  is  also  an  issue:  most  of 
the  experiments  are  conducted  in  laboratory-scale  combustors.  Nonlinear  processes,  especially  in  the  flame 
models,  are  not  completely  understood,  and  the  models  normally  include  a limited  representation  of  those 
processes.  In  the  case  of  laboratory  combustors,  it  is  possible  to  ignore  the  nonlinearities  and  design  a 
controller  on  the  system  linearized  around  the  ‘unstable’  equilibrium  point  that  the  control  aims  to  reach.  In 
the  case  of  classical  controllers,  they  scale  linearly  with  dimensions,  but,  at  this  point,  it  is  not  clear  how  the 
nonlinear  processes  in  the  combustor  scale:  they  could  become  a major  issue  in  large  combustors.  Moreover, 
a clearer  understanding  of  nonlinearities  is  required  for  the  design  of  nonlinear  controllers,  that  have  the 
potential  of  being  much  more  efficient  in  terms  of  the  ratio  of  the  energy  going  into  the  control  effort  versus 
the  total  energy  released  into  the  system.  Instabilities  in  combustion  chambers  typically  manifest  themselves 
as  pressure  oscillations  growing  to  limit  cycles.  Experiments  (Isella,  Seywert,  Culick  and  Zukoski  1997)  show 
that  in  some  combustors  there  is  a subcritical  bifurcation  leading  to  the  instability:  in  this  case  nonlinear 
control  might  be  more  appropriate,  and  results  from  linear  simulations  could  be  misleading,  especially  in 
the  vicinity  of  the  bifurcation  point.  For  the  case  of  heat  release  depending  on  velocity  fluctuation,  Wang 
(2000),  through  linear  analysis,  finds  conditions  on  the  heat  release  function  that  determine  if  the  bifurcation 
is  subcritical  or  supercritical.  In  particular,  he  has  shown  that,  if  the  heat  release  function  presents  a 
saturation  level,  the  amplitude  of  the  resulting  limit  cycle,  away  from  the  bifurcation  point,  is  independent 
of  the  actual  shape  of  the  heat  release  function.  This  result  plays  an  important  role  in  simplifying  the 
combustion  models  for  control  simulations.  If  the  conditions  for  supercriticality  are  satisfied,  using  the 
appropriate  value  for  the  saturation  of  the  heat  release  function  is  enough  to  capture  the  correct  dynamical 
behavior  of  the  combustor. 

While  several  of  the  studies  cited  include  sensor  noise  and/or  parameter  variations  as  part  of  the  system, 
none  of  them  makes  a clear  distinction  between  those  uncertainties,  the  intrinsic  noise  sources  of  the  system 
(additive  and  multiplicative),  and  the  unmodeled  dynamics.  In  the  following  we  will  describe  how  each  one  of 
these  effects  can  be  included  in  the  control  design  process.  Time-delay,  often  overlooked  in  the  consideration 
of  combustion  systems,  can  be  introduced  in  the  control  design  as  a further  uncertainty  by  using  the  same 
framework.  A different  method,  based  on  an  external  compensation  network,  is  also  presented. 

The  controller  is  designed  and  tested  on  the  model  of  a small  cylindrical  combustor  with  an  instability  in 
the  first  longitudinal  mode.  This  particular  model  has  been  widely  used  in  literature  as  a test  case  for  the 
analysis  (Culick  and  Yang  1992,  Yang,  Sinha  and  Fung  1992,  Haddad  et  al.  1997);  even  if  this  might  not  be 
realistic  in  the  specific  value  of  the  parameters,  the  model  presents  a dynamical  behavior  representative  of 
a combustion  chamber,  and  serves  well  as  a test  case.  The  methods  described  Me  obviously  not  limited  to 
this  specific  case,  and  can  be  applied  in  general. 

As  shown  below  the  noise  sources  arise  from  nonlinear  fluid  mechanics  and  as  such  form  an  integral  part  of  the 
system;  previous  studies  have  only  taken  into  account  external  noise  sources  (such  as  noisy  sensors/actuators). 
Recently  we  have  shown  how  these  intrinsic  noise  sources  affect  the  system  response  and  how  they  can  be 
used  to  identify  the  linear  parameters  of  a stable  system  (Seywert  and  Culick  2000). 
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2 Combustion  Instability  Model 


The  fluid  dynamical  equations  (continuity,  momentum  and  energy)  governing  the  flow  in  the  combustion 
chamber  can  be  combined  to  yield  a nonlinear  wave  equation  for  the  pressure  which  in  turn,  after  applying 
‘spatial  averaging’  leads  to  a system  of  coupled  oscillator  equations.  This  system  can  be  truncated  to  a finite 
number  of  modes  to  get  a low  order  representation  of  the  combustion  chamber.  This  procedure  has  been 
described  elsewhere  in  full  detail  (Culick  1976,  Culick  and  Yang  1992);  thus  only  a brief  summary  is  given 
here. 


We  begin  with  the  inhomogeneous  wave  equation  and  its  boundary  condition  for  the  pressure  fluctuation: 


= h 


1 d2p' 
a2  dt2 
n ■ Vp'  = 


(1) 


The  functions  h and  / contain  all  relevant  physical  processes  including  motion  of  the  boundary,  but  whereas 
the  contributions  from  gasdynamics  are  known  explicitly  every  other  process  included  requires  modeling,  a 
separate  matter.  Explicit  forms  for  h and  / are  given  in  (Culick  and  Yang  1992)  and  need  not  be  repeated 
here. 


Next  ‘spatial  averaging’  is  applied,  the  idea  being  that  any  unsteady  disturbance  in  a combustion  chamber 
can  be  synthesized  of  an  infinite  set  of  basis  functions,  chosen  normally  to  be  acoustic  modes.  For  a rocket 
motor  having  a chocked  exhaust  nozzle,  the  modes  are  those  for  a volume  enclosed  by  a rigid  boundary 
having  the  same  shape  as  the  internal  surface  of  the  motor,  but  with  no  combustion  and  flow. 

Hence  we  write  the  familiar  series  representation  of  the  pressure  field, 

OO 

p\x;t)=pJ2rij{t)il>j{x)  (2) 

i=i 

where  rjj  (t)  is  the  time-dependent  amplitude  of  the  jth  mode,  The  spatial  distribution  or  mode  shape  Y’j(x) 
is  calculated  as  the  solution  to 


= 0 
n • Vipj  — 0 


(3) 


and  kj  = aij  /a  is  the  wavenumber,  a is  the  mean  speed  of  sound  and  w j is  the  natural  frequency  of  the  mode. 


Substitution  of  (2)  in  the  left-hand  side  of  (la),  multiplication  of  the  equation  by  xj} „ and  integration  over  the 
volume  (i.e.  spatial  averaging  with  the  weighting  function  xjin)  leads  eventually  to  the  system  of  oscillator 
equations, 


cPtl n 
dt 2 


"b  ^nVn  ‘ — -Eri 


(4) 


where 


F„ 


xpnhdV  + 


and 


(5) 


The  central  idea  motivating  the  structure  of  the  analysis  is  that  combustion  instabilities  are  dominated 
by  acoustic  waves.  Hence  the  pressure  field  has  been  represented  by  the  expansion  (2)  with  p'  identified 
as  the  acoustic  pressure.  However,  following  a principle  discussed  by  Chu  and  Kovasznay  (1956),  small 
disturbances  are  in  general  made  up  of  three  kinds  of  waves:  acoustic,  entropy  and  vortieity  waves.  Unlike 
previous  analyses  in  which  only  the  organized  oscillatory  acoustic  field  was  accounted  for  all  three  wave  kinds 
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were  retained  in  the  analysis  given  by  Burnley  (1996)  and  Burnley  and  Culick  (1999).  The  extra  waves  give 
rise  to  stochastic  terms  in  the  equation  which  we  choose  to  retain  here. 

In  lowest  approximation  (small  amplitude,  uniform  mean  flow)  the  three  types  of  waves  propagate  indepen- 
dently and  hence  we  can  write: 


P 

0! 


Pa 

H'n 

4 

S'. + <$»  + *. 


(6) 


The  ‘forcing  function’  Fn  in  equation  (5)  is  a nonlinear  function  of  p'  and  u'.  In  previous  applications  of 
spatial  averaging  (Culick  and  Yang  (1992),  Culick  (1994)  and  references  cited  there)  the  velocity  fluctuation 
E is  related  to  p'  through  the  classical  linear  acoustics  equation: 


(7) 


With  entropy  and  vorticity  waves  present,  the  terms  ESI  and  v!3  make  new  contributions  to  F„.  Those 
contributions  can  be  written  explicitly  as  functions  of  En  and  Ea  (Burnley  and  Culick  1999)  but  since  no 
model  for  Ea  or  u'3  exists,  we  lump  them  together  as  stochastic  sources.  The  result  of  this  reasoning  is  the 
foEowing  system  of  acoustic  equations  that  are  the  system  used  in  previous  work  (CuHck  1995,  Culick  and 
Yang  1992,  Culick  1994)  with  the  additional  source  terms  representing  stochastic  or  noise  sources: 

CO  OO 

Vn  + ufcrin  = - + Fnity]  + 5^[&i(<)iji  + fn<(t)q»]  + S„(t)  (8) 

l=i  i=i 


Here  we  have  retained  only  the  Unear  terms  in  the  acoustic  amplitudes.  Since  we  are  going  to  design  a 
controUer  to  eliminate  the  pressure  oscillations  (i.e.  drive  all  acoustic  amphtudes  to  zero)  this  is  equivalent 
to  Unearizing  the  nonUnear  system  around  the  unstable  equilibrium  point.  Note  that  the  linear  terms 
include  all  linear  processes,  including  Unear  combustion  dynamics.  The  Unear  combustion  part  has  in  this 
formulation  been  lumped  together  with  the  linear  gasdynamics  into  the  coefficients  Eni  and  Dnv,  in  fact  it 
is  the  linear  combustion  that  makes  the  system  unstable  to  begin  with.  As  noted  by  (Wang  2000),  this  is 
a vaUd  approximation  as  long  as  the  system  presents  a subcritical  bifurcation  at  the  point  of  instability. 
Note  also  that  while  the  higher  order  acoustic  terms  have  been  neglected  the  nonUnear  dynamics  due  to  the 
vorticity  and  entropy  waves  are  included  in  the  noise  terms  f , and  S.  The  combustor  is  linearly  stable  if 

and  only  if  aU  modes  are  linearly  stable. 


For  control  appUcations  it  is  advantageous  to  reformulate  this  set  of  equations  in  state-space  form  with  state 
x = [qi  • - VnV i ■ ■ • Vn]t,  control  input  u,  and  output  y = pi  Ip.  The  foEowing  definitions  are  needed: 


0 1 

0 

-a2  - Enl  —D„i 

& = — 
p 

ip' actuator)  / 

C=[ip  i(  % sensor  )...1pN(x  sensor  )]  0) 


With  this  notation,  equation  (8)  - without  the  noise  terms  - becomes  equivalent  to 


x 

V 


Ax  + Bu 
Cx 


(10) 
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3 Control  Theory 

General  Considerations 

Even  before  the  development  of  models  including  combustor  dynamics  and  feedback  control,  experimental 
application  of  feedback  control  of  combustion  instabilities  was  successfully  tested  on  small  systems  (mainly 
using  loudspeakers  as  actuators).  Those  laboratory  demonstrations  report  examples  in  which  the  amplitudes 
of  limit  cycles  in  linearly  unstable  combustors  have  been  significantly  reduced,  sometimes  even  to  vanishingly 
small  values  (for  example:  Poinsot,  Bourienne,  Candel  and  Esposito  (1987),  Gulati  and  Mani  (1992)).  In 
most  cases,  the  ‘practical’  controller  was  a simple  proportional  feedback  or  a variation  of  a PID  controller. 
One  might  wonder  why  that  simple  approach  works  or,  conversely,  ask  why  we  need  more  sophisticated 
control  methods.  From  a general  viewpoint,  experiments  show  that  an  unstable  combustion  chamber  is 
a system  exhibiting  a linear  instability  (rapidly)  growing  to  a limit  cycle  (defined  by  the  non-linearities) 
that  typically  shows  a marked  predominant  frequency.  In  terms  of  dynamical  systems,  the  combustor  is 
characterized  by  two  unstable  complex-conjugate  poles  and  then  a series  of  stable  poles  with  relatively  large 
damping.  Provided  that  the  combustor  is  observable  and  controllable,  for  this  kind  of  system,  a proportional 
feedback  or  a PID  controller  can  be  successfully  tuned  to  obtain  a stable  feedback  loop  (Franklin,  Powell 
and  Emami-Naeini  (1995)).  Regarding  the  issue  of  controllability  (and  observability)  of  the  system,  for  the 
purpose  of  this  argument,  we  will  say  that  controllability  has  been  proved  in  practice  by  the  success  of  the 
experiments  cited.  A detailed  analysis  of  this  point  would  allow  optimization  of  the  position  of  actuators 
and  sensors,  but  that  is  out  of  the  scope  of  the  present  discussion.  The  need  for  more  sophisticated  control 
methods  derives  mainly  from  two  aspects:  first  one  might  want  to  impose  performance  specifications  on  the 
controller,  for  example  on  the  maximum  control  action,  or  on  the  noise  or  disturbance  rejection.  Second, 
combustion  systems  show  a high  degree  of  uncertainty  and  variability  (Lieuwen  and  Zinn  2000),  and  a 
controller  ‘tuned’  on  a particular  operating  point  does  not  guarantee  a reliable  performance.  Modern  control 
design  methods  allow  for  the  introduction  of  this  kind  of  consideration  during  the  synthesis  of  the  controller. 

All  the  considerations  above  and  the  design  method  presented  in  the  following  section  are  based  on  a linear 
model  of  the  combustor.  On  the  other  hand,  the  real  system  is  manifestly  nonlinear:  the  main  indication 
of  that  is  the  fact  that  the  pressure  oscillations  in  the  combustion  chamber  rapidly  reach  a limit  cycle.  A 
complete  understanding  of  the  dynamics  of  the  combustor  would  allow  tracing  the  source  of  the  nonlinear 
behavior  observed  in  the  experiments  (limit  cycles,  hysteresis,  as  in,  for  example,  Isella  et  al.  (1997)  and 
Lieuwen  and  Zinn  (2000))  to  its  origin:  nonlinear  gasdynamics  or  nonlinear  combustion.  In  that  case 
nonlinearities  in  the  system  could  be  exploited  by  an  ‘ad  hoc’  form  of  (nonlinear)  control  to  overcome  the 
main  limitations  of  linear  control:  requirement  of  a relatively  high  control  effort  and  actuation  frequency  at 
the  same  frequency  of  the  instability.  Since  such  a complete  model  is  not  available,  we  decided  to  limit  the 
analysis  to  the  linear  case.  Note  that  the  linear  model  of  the  combustion  chamber  presented  in  the  previous 
section  is  actually  a linearization  of  the  full  model  around  the  operating  point.  Since  the  main  purpose  here 
is  to  keep  the  system  ‘stable’,  i.e.  as  close  as  possible  to  the  linearized  equilibrium  point,  the  linear  model 
and  simulation  is  a valid  and  realistic  approximation  to  the  real  case,  provided  that  the  nonlinearities  do  not 
give  rise  to  a subcritical  bifurcation  (Wang  2000).  Note  that  nonlinearities  have  the  effect  of  limiting  the 
amplitude  of  the  oscillations:  hence  the  linear  model  is  in  this  sense  a ‘conservative’  approach  to  the  problem 
(for  example,  in  terms  of  required  control  action,  we  will  find  an  upper  limit).  In  short,  within  the  present 
approach,  nonlinearities  can  actually  be  neglected,  except  as  a formal  vehicle  for  rigorously  introducing  noise 
sources.  As  a consequence,  we  will  not  be  able  to  capture  the  effects  of  any  instability  mechanism  different 
from  the  linear  growth  and  phase  shifting  included  in  the  model  presented  above.  On  the  other  hand,  the 
present  approach  allows  for  a clear  distinction  of  the  effects  of  uncertainties,  intrinsic  noise  sources,  external 
noise  sources,  unmodeled  dynamics  and  time-delay. 
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Robustness 

The  controller  strategy  proposed  in  this  section  is  based  on  results  described  in  Chou,  Chen  and  Chao 
(1998)  and  Biswas  (1998).  It  allows  to  treat  parameter  uncertainty  (due  to  uncertainties  in  the  modeling  or 
system  identification  process),  multiplicative  noise  (in  this  case  intrinsic  to  the  system,  arising  from  vorticity 
and  entropy  waves)  and  residual  dynamics  (neglected  in  the  control  design  in  order  to  achieve  a low  order 
controller)  in  a unified  way. 

Including  parameter  uncertainty  in  the  controller  design 

Chou  et  al.  (1998)  considers  the  following  system: 

xc(t)  = [Ac  + AAe(f)]a:c(t)  -t-  AcrXr(t)  + [Bc  -)-  ABc(f)]u(f) 
xr(t)  = Arcxc(t)  + Arxr(t)  + Bru(t) 

V = [Cc  + A Cc(t)]xc  + CrxT  (11) 

In  this  description  the  system  has  been  split  into  two  parts:  the  controlled  dynamics  (state  xc)  which  will  be 
used  in  the  design  of  the  controller  and  the  residual  dynamics  (state  xr ) which  are  neglected  in  that  design. 
The  reasoning  behind  this  splitting  is  that  we  want  to  achieve  a controller  that  is  as  simple  as  possible.  This 
desire  leads  to  a need  for  a low  order  model  of  the  system;  the  controlled  dynamics  describe  that  low  order 
system  (which  at  the  minimum  needs  to  include  all  unstable  modes)  whereas  the  residual  modes  describe 
those  parts  of  the  original  system  that  the  designer  chooses  to  disregard  (the  higher  acoustic  modes  which 
are  strongly  attenuated  in  the  combustion  chamber). 

The  uncertain  parameters  of  the  controlled  system  are  assumed  to  be  bounded  (|  | denotes  the  modulus 
matrix): 

|A^c(i)|<QW,  |A£c(t)|  < and  |ACe(t)|  < Q(C)  (12) 

Q(A),  and  are  nonnegative  constant  matrices  that  describe  the  highly  structured  uncertainty  of 
the  parameters. 

Now  consider  a controller  based  on  a Kalman  filter  (used  to  reconstruct  the  state  from  the  measured  output, 
i.e.  the  pressure  signal  from  the  sensor)  with  estimator  gain  L and  feedback  gain  K: 

£c(*)  = Acxc{t)  + Bcu(t)  + L[y(t)  - CcS:x(t)} 

u(t ) = -Kxc(t)  (13) 

Then  the  controlled  closed-loop  system  is  described  by  Hc  and  its  interaction  with  the  residual  system  is 
governed  by  Hcr  and  Hrc  (note  that  Hcr  and  Hrc  have  been  changed  to  allow  for  the  presence  of  Acr  and 
Arc  when  compared  to  Chou  et  al.  (1998)). 

] Kr,,\K.-B,K  - B,K)  (14) 

Define  the  matrices  Gc  and  Gr  as  follows: 

Gr  > [ sup  I3t^,(jw)|  ] and  Gc  = [ sup  |p^(jw)|  ] (15) 

Where  u > 0 and  and  are  the  ikth  elements  of  (ju/  — Ar)~1 , respectively  (jail 
the  uncertainty  matrix  U as 

q^  + qW\k\  qW\k\  ■ 

“ L Q(A)  + Q(A)\K\  + \L\Q(C ) QW\K\  _ 


— Hc)  1 . Also  define 
(16) 


According  to  results  given  in  Chou  et  al.  (1998)  (with  the  trivial  extension  to  include  Acr  and  Arc),  the 
closed-loop  system  will  be  stable  if  the  matrices  Hc  and  Ar  are  stable  matrices  and  the  following  inequality 
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Combining  both  approaches  and  designing  a controller 

We  can  combine  both  approaches  (i.e  account  for  bounded  parameter  uncertainty  and  random  Gaussian 
noise  perturbations)  by  including  the  H-terms  in  equation  (20)  in  the  uncertainty  of  the  closed-loop  matrix 
Hc  in  equation  (3).  Thus  by  redefining  the  uncertainty  matrix  U, 

Q<*>+Q<B>|K|+2<AA>+3(MB>|/r|+(H<XB>|ir|)r+|K|T=<BB>|K|+|z,|Ts<cc>|i|  Q<B>|*r|+5<AB>|iC|+|K|T=(BB>|*:| 

— Q<A'+<3<A>|iC|+|i-|Q(c>+{H(AB>|K-|)T+|JC|'rS(Be>|/C|  Q<b>|K-|+|*-|tS<bb>|K| 
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and  fulfilling  equation  (17),  robust  stability  can  still  be  guaranteed. 

Up  to  this  point  we  have  derived  conditions  that  L and  K need  to  fulfill  without  specifying  how  to  design 
them.  Note  that  the  controller  is  to  be  designed  for  the  nominal,  undisturbed  plant;  we  can  then  check 
equation  (17)  to  make  sure  it  is  also  effective  on  the  perturbed  system  - or  we  can  use  that  equation  to  see 
how  much  noise  or  parameter  variation  our  controller  can  handle. 


The  design  of  the  controller  depends  of  course  on  the  performance  we  wish  to  achieve.  Any  method  can  be 
used:  Biswas  (1998)  uses  a pole-placement  technique,  here  we  follow  Chou  et  al.  (1998)  in  using  a standard 
LQG  method. 

LQG  is  advantageous  because  it  allows  for  the  inclusion  of  additive  system  noise  (£(x>)  as  well  as  sensor  noise 
(£<»)).  Thus  the  complete  system  we  try  to  control  is  given  by: 

ic(t)  = [Ac  + A Ae(t)  + tlA)(t)}xc{t)  + Acrxr(t)  + [Be  + ABc(t)  + S{B)(t)}u(t) 

Xr(t)  = ArcXc(t)  + ArXr(t)  + Bru(t)  + ^ 

y = [Cc  + ACc(t)  + ^C)(t)}xc  + Crxr+  (21) 


In  using  this  method  we  minimize  the  performance  index  J given  by 

roo 

J-  [xJ(t)Qxc(t)  + pcuT(t)Ru(t)]  dt 

Jo 


(22) 


In  this  case  pc  is  the  design  parameter  that  will  be  changed  to  fulfill  condition  (17).  Since  we  want  to  reduce 
the  pressure  oscillations  at  a specific  location  xp,  it  is  natural  to  choose  Q such  that  Qxc  represents  our 
best  estimate  of  p'  at  this  point 


Q=P2- 


V’i(zp)  o o 

o 0 

0 0 t^(Xp)  . 


and  R = 1 


(23) 


Time  Delay 

Time  delays  often  arise  in  combustion  systems:  for  example,  even  when  no  control  is  present,  there  is  delay 
between  injection  of  the  fuel  mixture  and  fully  developed  combustion  for  the  case  of  liquid  or  gas  combustors. 
When  feedback  control  is  present,  there  are  delays  intrinsic  to  the  controller  due  to  finite  rates  of  actuators 
and  sensors,  time  spent  for  signal  acquisition  and  processing,  and  clock  time  in  case  a digital  computer  is 
used.  Even  for  the  typical  laboratory-scale  combustor,  when  a loudspeaker  is  used  as  an  actuator,  time 
delays  might  play  an  important  role:  suppose  the  first  unstable  acoustic  mode  has  a frequency  of  1 kHz, 
then  a typical  reaction  time  for  the  controller  (if  we  consider  linear  approach,  the  bandwidth  of  the  controller 
should  at  least  match  the  typical  frequencies  of  the  instability  in  the  plant)  is  of  the  order  of  1ms.  Modem 
electronic  equipment  can  certainly  process  the  required  computation  for  determining  the  control  input  much 
faster  than  that;  the  bottleneck  for  this  case  is  the  time  it  takes  for  the  pressure  input  (from  the  loudspeaker) 
to  influence  the  chamber  acoustic  response.  This  time,  for  a 50cm  chamber,  is  of  the  order  of  l-2ms,  just 
the  same  order  of  the  instability.  In  the  case  of  industrial  scale  combustors,  or  when  using  secondary  fuel 
injection  as  control  actuation,  the  necessity  of  considering  time  delays  becomes  even  more  compelling,  since 
in  these  cases  the  time  delay  can  easily  be  larger  than  the  characteristic  timescale  of  the  instability. 

Time  delays  always  reduce  the  stability  of  a system  (Franklin  et  al.  1995),  hence  it  is  very  important  to 
take  them  into  consideration  when  simulating  a realistic  combustor  and  when  designing  a suitable  controller. 
Regarding  the  controller  design  phase,  three  general  approaches  are  possible. 
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• Classical  Control.  If  we  look  at  the  transfer  function  of  the  system,  and  indicate  with  r the  time  delay, 
the  problem  with  time  delay  is  reduced  to  a conventional  one  by  expressing  the  non  rational  function 
eT’  in  terms  of  a rational  function.  Note  that  the  function  eTS  is  analytic  (for  finite  values  of  s),  so 
approximation  with  a rational  function  is  allowed.  A typical  approach  is  to  use  a Pade  approximant, 
based  on  a McLauren  series  expansion  of  the  exponential  function.  The  value  of  the  method  is  limited 
by  two  factors.  First,  the  rationed  approximation  of  the  delay  rapidly  increases  the  effective  order  of 
the  plant,  making  the  control  design  problematic.  Second,  large  values  of  time  delay  will  decrease  the 
available  phase  margin  to  the  point  where  it  is  no  longer  possible  to  design  a stabilizing  controller. 
Also,  (Wang  2000)  shows  that  a low  order  polynomial  approximation  of  the  time  delay  is  not  enough 
to  have  a satisfactory  model  in  terms  of  dynamical  behavior  of  the  original  system.  In  the  present 
paper  we  do  not  take  this  method  into  consideration,  since  we  focus  our  attention  on  control  design 
methods  capable  of  incorporating  robustness  requirements. 

• Modem  Control.  In  this  case,  time  delay  can  be  viewed  as  an  uncertainty  in  the  system  and  incorporated 
in  the  design  as  a perturbation  to  the  original  plant.  More  details  are  given  in  the  following  section. 

• Delay  Compensation.  This  category  includes  all  the  other  methods  used  to  compensate  for  time  delay. 
An  important  group  includes  compensation  networks  that  bring  the  delay  ‘out  of  the  loop’,  and  hence 
allow  to  design  the  controller  using  conventional  methods  applied  to  the  plant  without  time  delay.  A 
typical  example  is  the  Smith  Regulator.  A caveat  here  is  constituted  by  the  fact  that  most  networks 
based  on  linear  elements  generally  do  not  modify  the  eigenvalues  of  the  original  plant,  so  they  only 
apply  to  stable  (or  marginally  stable)  plants.  On  the  other  hand,  by  using  these  methods,  arbitrarily 
large  time  delays  can  be  accounted  for  without  loss  of  stability  margin. 


Modern  Control  Approach 

Time  delay  can  be  incorporated  in  the  design  of  a controller  by  considering  the  time  delay  as  a multiplicative 
perturbation  to  the  plant.  Let  P( s)  = P(s)e~T>  be  the  perturbed  plant.  The  perturbed  plant  can  be  included 
in  the  set 

{(1  + A piantWunc)P  :||  Aplant  |U<  1}  (24) 

with  the  choice  of  an  appropriate  weight  function  Wurlc 

\e-Tju-l\<\Wunc{jw)\  Vu>,r  (25) 

The  design  of  the  controller  then  proceeds  in  the  same  way  as  before.  Note  that,  if  the  time  delay  is  ‘large’, 
condition  (25)  typically  imposes  a significant  limitation  on  the  controller;  in  general  performance  is  degraded 
and,  if  r is  large  enough,  it  might  be  impossible  to  design  a stabilizing  controller  for  the  delayed  system.  On 
the  other  hand,  when  a solution  exists,  stability  and  performance  are  guaranteed  according  to  the  design. 
Uncertainty  in  the  numerical  value  of  the  time  delay,  as  it  is  typical  in  combustion  systems,  is  automatically 
taken  into  consideration  by  the  design  method.  Application  of  this  approach  is  included  in  the  example 
presented  later. 


Delay  Compensation  Approach 

In  this  section  we  will  examine  a method  based  on  predictive  control:  the  time  delay  is  compensated  by  a 
predictor  that  acts  on  the  measured  or  estimated  state  and  feeds  the  controller  with  the  appropriate  signal 
to  perform  the  feedback  action  at  the  compensated  time.  The  control  system  consists  of  a predictor  and  a 
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controller;  the  closed  loop  equations  can  be  written  as 

x(t)  = Ax(t)  — Bu{t  — h) 
y(t)  = Cx(t) 

r° 

p(i)  = eArx(t)  + e ATBu(t  + t)dt 
u{t)  = Kp[t)  + i{t) 

where  i(t)  is  the  external  input  to  the  system  and  might  not  be  present.  The  predictor  written  as  in  (27)  is 
simply  derived  by  integrating  (26)  from  the  current  time  t to  the  time  t + r.  A change  of  variables  produces 
the  form  (27),  which  contains  information  only  up  to  the  current  time  and  consequently  it  can  be  physically 
implemented. 

Furukawa  and  Shimemura  (1983)  prove  (based  on  direct  computation  of  the  closed-loop  characteristic  equa- 
tion) that  if  the  pair  ( A,B ) is  controllable,  then  the  predictor  (27)  and  the  controller  (28)  yield  a finite 
spectrum  of  the  closed-loop  system,  located  at  arbitrarily  preassigned  points  in  the  complex  plane. 

Note  that  the  predictor  (27)  contains  an  integral  term  up  to  the  current  time.  It  is  impossible  to  integrate 
up  to  current  time  without  solving  an  integral  equation,  or  iterating  on  the  solution,  but  it  can  be  shown 
(Manitius  and  Olbrot  1979)  that  the  limits  in  the  integral  term  appearing  in  (27)  can  be  substituted  by 
— r — e and  — e if  e is  sufficiently  small.  Let  us  now  consider  the  following  scheme: 


(26) 

(27) 

(28) 


i(t)  = Ax(t)  — Bu(t  — h) 

(29) 

z{t)  = Az(t)  — Bu(t) 

z(t)  = z(t)  — eArz(t  — t) 

(30) 

p(t)  = eArx(t)  + z(t ) 

it(f)  = Kp{t) 

(31) 

A simple  substitution  shows  that  computing  z(t)  from  the  equations  above  results  in  the  evaluation  of  the 
correct  predictor  term  (27).  The  scheme  (29)-(30)  can  be  physically  implemented  in  Simulink  by  using 
the  network  connection  presented  in  figure  1,  where  symbols  refer  to  the  letters  used  in  (29)-(31)  and  p 
is  the  signal  sent  to  the  controller.  Since  the  time  delay  is  compensated  in  this  secondary  predictor-loop, 
application  of  this  scheme  to  the  plant  presented  above  allows  to  include  significant  time  delays  in  the  system 
without  compromising  performance  in  the  design  of  the  controller.  On  the  other  hand,  the  use  of  a second 
loop  reduces  the  robustness  of  the  system  to  uncertainties  in  the  value  of  the  parameters. 


O 
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Figure  1:  SIMULINK  realization  of  the  predictor  block,  equations  (30) 
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Ail 

71  = 1 

n — 2 

n = 3 

n = 4 

Enl 

71=1 

n = 2 

n = 3 

n = 4 

1=1 

-.01 

.007 

-.001 

.007 

1=1 

-.005 

-.005 

.0025 

.0016 

1 = 2 

.01 

.1 

.007 

-.001 

1 = 2 

-.0025 

-.015 

.01 

.01 

1 = 3 

-.01 

.01 

.75 

.008 

1 = 3 

-.005 

.0 

-.02 

.02 

1 = 4 

.02 

-.005 

.01 

1.5 

1 = 4 

.01 

.02 

.02 

-.025 

Table  1:  Combustion  chamber  parameters 


4 Example 


In  the  following  we  will  briefly  demonstrate  the  described  design  procedures  on  a particular  example.  The 
example  has  been  chosen  solely  because  it  has  been  used  previously  in  the  literature.  The  methods  used 
are  obviously  much  more  general  and  can  be  applied  to  tiny  linearized  (combustion)  system  - as  mentioned 
previously,  all  processes,  including  combustion,  have  been  linearized  and  are  embedded  in  the  model  para- 
meters. 


The  numerical  example  used  is  the  same  one  as  given  in  Haddad  et  al.  (1997)  and  used  previously  by  Culick 
and  Yang  (1992)  and  Yang  et  al.  (1992).  The  combustion  chamber  is  assumed  to  be  cylindrical  of  length  L 
and  only  longitudinal  modes  are  considered.  The  chamber  is  closed  on  the  upstream  end  and  has  a nozzle 
at  the  downstream  end  which  acts  as  an  acoustically  closed  boundary.  The  sensor  detecting  the  instability 
is  a microphone  located  at  x,  whereas  the  actuator  used  to  control  the  pressure  oscillations  is  a loudspeaker 
located  at  xa . The  internal  dynamics  of  the  loudspeaker  are  modeled  as  a second  order  system  and  included 
in  the  state-space  formulation,  i.e.  they  form  an  integral  part  of  Ac  in  equation  (21)  through  augmentation 
of  the  state  xc.  In  this  way  the  actuator  dynamics  are  accounted  for  in  the  design  of  the  controller.  In  this 
example  the  actuator  (and  sensor)  are  treated  as  perfectly  known  systems;  within  the  framework  described 
here  it  is  straightforward  to  include  uncertainties  or  noise  in  those  systems  too.  Note  that  the  model  of  the 
actuator  as  a ‘loudspeaker’,  a second  order  system  with  a high  gain,  is  actually  more  general  than  it  seems. 
If  we  wanted  to  represent  an  injector,  the  same  model  would  still  apply,  with  a longer  time  delay,  and  some 
difference  in  the  numerical  value  of  the  parameters,  but  substantially  the  same  structure,  i.e.  second  order 
dynamics  and  very  high  gain;  see  Neumeier,  Nabi,  Arbel,  Vertzberger  and  Zinn  (1997)  for  an  example. 

The  linearized  model  of  the  combustion  chamber  is  characterized  through  the  parameters  Eni  and  D,d  which 
are  given  (after  non-dimensionalization  of  time  by  ^)  in  table  1 for  the  first  4 longitudinal  modes  - the 
parameters  are  originally  taken  from  Yang  et  al.  (1992)  and  are  typical  for  solid  propellant  rockets. 


In  this  design  study  we  only  consider  uncertainty  in  the  parameters  E,d  and  D„i  and  noise  in  those  terms 
given  in  equation  (8).  In  other  words  there  is  no  noise  or  uncertainty  in  the  sensing  or  actuation  process 
and  the  stochastic  sources  act  (as  described  by  the  model)  as  random  perturbations  of  Eni  and  D„i.  Thus 
we  put  (keeping  the  notation  of  the  previous  section) 

Q(B)  = 0 Q(c)  = 0 E^AB)  = 0 xSBB)  - 0 E(CC)  = 0 (32) 


Furthermore  we  assume  that  all  stochastic  sources  (due  to  the  vorticity  and  entropy  modes  in  the  chamber) 
are  uncorrelated  and  have  the  same  variance  a2.  Fin  silly  we  assume  that  the  parameter  uncertainties  can  all 
be  described  by  a single  variable  e.  These  assumptions  are  only  made  to  simplify  the  expressions  as  we  can 
now  write: 


QW  = 6 • 


0 0 
I Ail  I Aii| 


and  =N-cr2  ■ 1 


(33) 


A controller  was  designed  using  an  LQG  technique  by  taking  only  1 mode  (the  unstable  first  mode)  into 
consideration.  It  is  assumed  that  the  complete  system  is  given  by  4 modes  and  thus  the  remaining  3 modes 
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Figure  2:  Left:  Nominal  system  (only  1 mode)  response  to  control.  Right:  Response  of  complete  system  (all 
4 modes)  to  the  same  controller 


are  considered  to  give  the  residual  dynamics.  Basing  the  controller  on  a minimal  set  of  modes  is  desired 
as  it  reduces  reduces  the  order  of  the  controller  and  thus  allows  for  easier  implementation.  Figure  2 shows 
the  response  of  the  reduced  and  complete  system  to  the  controller  (turned  on  at  t = 100).  As  expected  the 
presence  of  the  extra  modes  (not  considered  in  the  design)  in  the  full  system  reduces  the  performance  of  the 
controller  (slower  decay). 

Figure  3 shows  the  guaranteed  stability  limits  (in  terms  of  e and  a)  of  the  controller.  The  solid  line  is  the 
limit  predicted  by  equation  (17)  for  the  truncated  system  where  only  1 mode  is  used  in  the  simulation. 
The  other  lines  describe  the  stability  region  as  more  modes  are  added  in  the  simulation,  i.e.  as  the  system 
approaches  the  ‘complete’  system.  The  region  shrinks  as  additional  modes  are  introduced  into  the  simulation 
while  the  same  controller  (based  on  1 mode)  is  retained.  In  this  extreme  case  (where  we  considered  only  one 
mode  to  base  the  controller  on)  the  changes  are  drastic,  but  as  the  neglected  modes  become  more  heavily 
damped  their  influence  grows  smaller  (as  can  be  seen  by  the  lines  moving  closer  together). 

Figure  4 illustrates  the  effect  the  (stable)  residual  system  can  have  in  the  presence  of  noise.  The  same 
controller  (based  only  on  the  unstable  first  mode)  is  used  in  all  three  cases.  In  case  A (low  noise)  the 
controller  is  able  to  stabilize  the  full  system  (all  4 modes  included  in  the  simulation);  however  at  a higher 
noise  level  (case  B)  the  pressure  oscillations  do  not  decay  to  zero.  Note  that  this  noise  level  is  well  within 
the  stability  limits  as  predicted  with  the  truncated  (1  mode)  system  and  thus  underlines  the  importance  of 
the  neglected  dynamics.  If  the  simulation  is  performed  with  the  reduced  system,  the  instability  does  indeed 
decay  as  anticipated  (case  C). 


Example  with  Time  Delay  Compensation 

The  delay  compensation  approach  allows  separating  the  control  problem  from  compensation  of  the  time 
delay.  Figure  5 presents  the  results  of  the  application  of  the  method  to  our  model  combustor.  The  non- 
dimensional  time  delay  is  chosen  to  be  r = 10,  which  corresponds  to  a delay  of  about  10ms,  i.e.  5 periods 
of  an  oscillation  at  500 Hz  and  constitutes  a reasonable  upper  limit  to  the  delay  that  can  be  expected  in  a 
real  combustor  controlled  by  modulating  the  injection  of  a secondary  fuel.  Note  how  the  predictor  works: 
the  controller  (control  action  is  plotted  in  the  bottom  half  of  figure  5)  starts  sending  commands  immediately 
when  it  is  activated.  The  control  is  computed  on  a prediction  of  the  future  state  of  the  system,  i.e.  the  state 
of  the  system  when  the  control  signal  will  effectively  reach  the  plant.  The  system  response,  plotted  in  the 


B,C 


unstable 


A-> 


1 mode  controller 


Figure  3:  Guaranteed  stability  limits  for  dosed-loop  system  in  a — e plane  (noise  variance  and  uncertainty) 
depending  on  number  of  residual  modes  considered 


Figure  4:  From  left  to  right:  time  simulations  of  cases  A,B  (complete  system),  and  C (reduced  system)  as 
defined  in  figure  3. 


top  half  of  5,  shows  that  the  system  effectively  starts  reacting  to  the  control  at  a non-dimensional  time  of 
40  when  the  controller  is  put  online  at  a non-dimensional  time  of  30. 


5 Concluding  Remarks 


In  this  paper  we  have  shown  how  uncertainties,  noise,  unmodeled  dynamics  and  time-delay  can  be  included 
in  the  controller  design  for  combustion  instabilities. 

A clear  distinction  has  been  made  between  the  uncertainty  and  the  noise.  This  is  necessary  as  the  parameter 
uncertainty  can  be  bounded;  e.g.  in  practical  applications  we  might  know  that  in  the  operating  range  of 
interest  the  various  parameters  are  located  within  certain  numerical  bands.  In  contrast,  true  noise  sources 
can  in  general  not  be  bounded,  and  thus  do  not  fit  in  the  common  control  frameworks;  they  are  characterized 
by  their  mean  values,  which  we  include  in  the  system  parameters  and  by  their  variation. 
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Figure  5:  Control  with  time-delay,  (r  = 10).  Top  half:  system  response.  Bottom  half:  control  action. 


Explicit  consideration  of  the  neglected  modes  allows  studying  their  influence  on  the  controller  robustness. 
This  is  especially  important  since  in  most  experimental  implementations  to  date  the  controller  has  been 
designed  by  taking  only  the  unstable  mode(s)  into  account.  In  the  example  given  here  only  the  first  mode 
is  unstable  and  it  is  in  fact  possible  to  stabilize  the  system  by  solely  controlling  this  one  mode.  Note  that 
the  controller  is  designed  to  accommodate  large  uncertainties  (or  noise)  since  we  anticipate  that  the  residual 
modes  will  affect  the  dominant  first  mode.  This  is  the  way  unmodeled  dynamics  are  traditionally  handled: 
by  including  them  in  the  uncertainties  of  the  system.  The  framework  presented  here  shows  how  much  of  that 
uncertainty  can  be  attributed  to  the  neglected  modes.  In  the  example  given,  the  damping  of  the  ignored 
modes  (notably  the  second  one)  is  rather  small  and  thus  we  see  that  the  actual  parameter  uncertainty  e (or 
noise  intensity  cr)  that  the  controller  can  tolerate  declines  dramatically  as  additional  modes  are  considered. 
Therefore  we  conclude  that  the  residual  dynamics  dominate  the  uncertainty  unless  the  neglected  modes  are 
highly  damped. 

Inclusion  of  a time-delay  in  the  modern-design  framework  as  an  uncertainty  is  adequate  when  the  time-delay 
is  of  the  same  order  of  the  characteristic  time  of  the  instability,  defined  as  the  inverse  of  the  frequency  of 
the  unstable  mode.  Cases  with  longer  time  delays,  as  it  might  be  the  case  in  full-scale  combustors,  can 
be  treated  by  adding  a second  loop  to  compensate  for  the  delay:  simulation  shows  very  good  performance, 
but  issues  about  robustness  to  uncertainty  and  perturbations  need  to  be  addressed  carefully.  An  adaptive 
observer  might  be  needed  for  application  to  real  systems.  Future  work  in  this  area  should  include  testing 
of  the  concepts  on  an  experimental  combustor;  and  system  identification  to  define  better  models  of  real 
actuators,  in  particular  injectors  and  fuel  flow  modulators.  More  analysis  is  also  needed  to  characterize  (and 
eventually  take  advantage  of)  nonlinearities  naturally  present  in  combustion  chambers. 
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PAPER  -27,  F.  Culick 

Question  (F.  Vuillot,  France) 

To  analyze  numerical  simulation  results  in  terms  of  control  efficiency,  it  would  be  useful 
to  apply  the  flow  decomposition  (into  acoustic-vorticity-entropy  fields)  to  the  energy.  Do 
you  know  of  work  in  that  direction? 

Reply 

Not  specifically  for  the  situation  I discussed  here.  However,  I believe  that  the  basis  exists 
for  doing  so.  The  equations  exist,  of  course,  and  I suggest  that  work  by  Morfey  many 
years  ago  in  acoustic  energy  in  moving  media  would  be  a good  place  to  start.  Flancho  (?) 
later  used  Morfey’s  ideas  in  his  work  on  nonlinear  acoustics;  his  approach  and  results 
may  well  be  helpful. 
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